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Abstract 

We present two models for electron-driven uphill proton transport across lipid membranes, with 
the electron energy converted to the proton gradient via the electrostatic interaction. In the first 
model, associated with the cytochrome c oxidase complex in the inner mitochondria membranes, 
the electrostatic coupling to the site occupied by an electron lowers the energy level of the proton- 
binding site, making the proton transfer possible. In the second model, roughly describing the redox 
loop in a nitrate respiration of E. coli bacteria, an electron displaces a proton from the negative 
side of the membrane to a shuttle, which subsequently diffuses across the membrane and unloads 
the proton to its positive side. We show that both models can be described by the same approach, 
which can be significantly simplified if the system is separated into several clusters, with strong 
Coulomb interaction inside each cluster and weak transfer couplings between them. We derive 
and solve the equations of motion for the electron and proton creation/annihilation operators, 
taking into account the appropriate Coulomb terms, tunnel couplings, and the interaction with the 
environment. For the second model, these equations of motion are solved jointly with a Langevin- 
type equation for the shuttle position. We obtain expressions for the electron and proton currents 
and determine their dependence on the electron and proton voltage build-ups, on-site charging 
energies, reorganization energies, temperature, and other system parameters. We show that the 
quantum yield in our models can be up to 100% and the power-conversion efficiency can reach 
35%. 

PACS numbers: 82.39.Jn, 87.16.A-, 73.63.-b 
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I. INTRODUCTION 



Every living organism obtains the energy needed for its survival from the outside world. 
This energy can be in the form of sunlight or food; but in both cases it is unstable and cannot 
be utilized directly, so several energy-conversion steps are necessary. One of the most widely 
used intermediate forms for energy storage is the electrochemical proton gradient across lipid 
membranes, such as the inner mitochondrial membranes or plasma membranes in bacteria. 
To achieve and maintain this proton gradient, nature employs several different types of 
electron- or light-driven systems, where the energy of high-energetic electrons or absorbed 
photons is used for the energetically-uphill proton transfer from the negative (JV) to the 
positive (P) sides of the membrane. 

Here we discuss two mechanisms of energy conversion from the highly unstable electronic 
form of energy to the proton gradient, namely, proton pumps and redox loops [l, 2\. Both 
mechanisms rely on the electrostatic interaction between electrons and protons, although the 
specific details of the proton pumps and the redox loops look very different. For example, in 
a proton pump, such as cytochrome c oxidase, electrons move mainly along the membrane, 
whereas protons move across the membrane, which results in an accumulation of the positive 
charge on the P-side and in the generation of a proton- motive force (PMF) 3j-|6|. In the 
redox-loop mechanism of PMF generation, taking place in the nitrate respiratory chain 
of E. coli bacterium, the neutral shuttle, carrying both protons and electrons, crosses the 
membrane. Here, the charge accumulation occurs when electrons cross the membrane, just 
before embarking on the shuttle, and right after unloading from the shuttle Q-Q. It should 
be noted that the proton pump operating in the cytochrome c oxidase has no essential 
mechanically-moving parts, whereas the redox-loop mechanism is impossible without the 
molecular shuttle diffusing between the negative and the positive sides of the lipid membrane. 

In general, the treatment of the electron and proton transfer events is extremely difficult 
because the total number of the occupation states increases exponentially with the number 
of the electron- and proton-binding sites, when all of them are electrostatically coupled. In 
the present work, however, we show that both above-mentioned mechanisms of the trans- 
membrane proton translocation can be described with a similar mathematical model, taking 
into account the Coulomb interaction between one electron- and one proton-binding sites 
only, and neglecting electrostatic couplings to other sites. It is necessary to have at least 
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three redox sites and three proton-binding sites in order to obtain a proton pumping effect 
and suppress a reverse flow of protons from the P-side to the iV-side of the membrane. In 
the absence of strong Coulomb interaction between all sites, there is no need to introduce a 
electron and proton occupation states (as was done in our previous works, 
), which grows exponentially with the number of sites. Instead, we now di- 
vide the whole system into clusters of strongly coupled sites. These clusters are described by 
their own set of occupation states, and the total number of the states in the system is equal 
to the sum (not the product!) of the states in the clusters. The clusters are weakly coupled 
by electron tunneling terms and by proton transfer amplitudes, so that transitions between 
the clusters can be considered within perturbation theory. While in this work we present 
quite simple models, similar approaches can be applied to much more complicated biological 
systems, such as Photosystem II and the whole respiratory chain in the inner mitochondrial 
membrane j^j. 

The quantum yield for the two models analyzed in this paper can be about 1. Why such 
a high quantum yield? This can be explained from the fact that, in order to be transferred 
through the system, an electron needs to loose its energy. This cannot be done via the 
environment because the reorganization energy is not large enough. Consequently, electron 
transport occurs with the assistance of protons gaining this energy and being transferred to 
the positive side of the membrane. Thus, the transfer of a single electron is accompanied by 
the transfer of a single proton and the corresponding currents are equal, which results in an 
almost perfect quantum yield. 

II. MODEL 

We consider a physical model describing an electron-coupled translocation of protons 
from the negative (N) to the positive (P) side of a membrane. The model consists of an 
interaction site, Q = {Q e , Q p }, containing a single electron level with energy Eq and a single 
proton energy level characterized by the energy Eq. We also introduce two electron sites, 
L and R, coupled to the electron site Q e , and two proton sites, A and B, coupled to the 
proton site Q p (Fig. 1). The electron site L is coupled to the electron source S, and the site 
R is connected to the electron drain D. The proton site A is coupled to the proton reservoir 
N (the negative side of the membrane), and the site B is coupled to the positive side of the 
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membrane (proton reservoir P). 



A. Hamiltonian 

The Coulomb interaction between an electron and a proton, both located on the central 
site Q, is described by the energy uq, so that the Hamiltonian of the site Q has the form 

H Q = SqUq + EqNq - UqUqNq, (1) 

where uq = clqCiq is the electron population of the site Q, and Nq = 6q6q is the proton 
population of this site. Electrons are described by the Fermi-operators a a , and protons are 
characterized by the Fermi-operators b a with a = L, Q e , R and a = A, Q p , C, and with the 
corresponding populations n a = a\a a , N a = blp a . 

The contribution of the electron sites L, R and the proton sites A, B to the total Hamil- 
tonian of the system is described by the term 

H = e L n L + e R n R + E A N A + E B N B , (2) 

where El,£r are the energy levels of the electron sites L and R, and E A , E B are the energies 
of the proton-binding sites A and B. 

The strongly-interacting electron and proton sites Q e and Q p form a single (interaction) 
cluster, whereas the sites L, R and A, B separately form other four (peripheral) clusters. 
The cluster Q can be characterized by the vacuum (empty) state and by three additional 
occupation states, or, equivalently, by the average electron and proton populations, (uq) 
and (Nq), complemented by the correlation function, K = (uqNq). The other electron and 
proton clusters are described by the corresponding average occupations, (ni), {ur) and 
(N A ), (N B )- For six electron and proton-binding sites we should have 2 6 = 64 occupation 
states. However, with the cluster approach, the system can be completely described by 
only seven functions: (nq), (Nq), K (for the interaction cluster), and (ul), (hr), (Na), (Nb) 
(for the peripheral clusters). Previously, we applied a similar approach to analyze quantum 
transport problems in nanomechanical systems 14^ . 
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1. Electron and proton transitions 



The electron tunneling Hamiltonian between the site Q and the sites L and R is given by 

He — —Al a L a,Q — Ar a R <jq + H.c., (3) 

whereas the A-Q and B-Q proton transitions are described by the term 

H p = -A A b A b Q - A B b B b Q + H.c. (4) 

Here A L ,A R are the electron tunneling coefficients, and Aa, A B are the proton transfer 
amplitudes. In the case of a movable interaction site, e.g., when the electron and proton 
sites Q are located on the shuttle (quinone/quinol), the amplitudes A L ,A R and A A ,A B 
depend on the position x of the shuttle. 

The S-lead serves as a source of electrons, and the D-lead works as an electron drain. 
The coupling to these leads is characterized by the Hamiltonian 

H L r = - ^tks 4s a L ~ ^2 tkD c kD a R + H - c - ( 5 ) 

k k 

The proton transitions between the iV-side of the membrane and the site A and between 
the P-side of the membrane and the site B are described by the Hamiltonian 

H AB = -J2 t i n d U b A~Y. T i p d \p b ° + H - c - ( 6 ) 

Here Cks, CkD are Fermi operators of the electron reservoirs S and D, and d q ^, d q p are the 
Fermi operators of protons in the reservoirs iV and P. The electron reservoirs S and D have 
the Hamiltonian 

HsD = ^2( £ksC kS C kS + £kDcl D C kD ), (7) 
k 

and are characterized by the Fermi distributions fs(£ks), /d( £ *;d) with the corresponding 
electrochemical potentials /is and ^d- For the proton reservoirs iV and P we have the 
Hamiltonian 

Hnp = y^^EgN d qN d qN + E qP d qP d qP ), (8) 
q 

with the Fermi distributions Fjv(-E 9 iv) and Fp(E q p) and the proton electrochemical poten- 
tials /in and \ip. 
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2. Environment 



The interaction of the electron-proton system with the protein environment, which is 
described as a sum of independent oscillators (l5j, is characterized by the Hamiltonian 



,2 



H cnv - ^ + ~^T~ ( X i ~ X 3° n ° ~ X 0S X] C kS C kS - XjD ^ C kD C kD- 



2m j 

3 3 



^ Xj a N a - X jN d\ N d qN - X jP ^2 d\pdqP j 

a 9 q / 



(9) 



where n CT = aJ.a CT is the population of the electron site a [a = L,Q,R), N a = b\jb a is 
the population of the proton site a (a = A,Q,B). The constants Xj a ,Xj$,XjD determine 
the electron coupling to the environment, and the parameters Xj a , Xj^, Xjp describe the 
proton-environment interaction. 
With the unitary transformation, 



U = exp 



J2 P j(j2 X i° U ° + X i S + XjD C kD C kD + 

^ cr k k 

Xj a N a + Xj N cfi qN d q N + Xjp dqpd q p 



i 

3 



and 



with the phases 



& = J2Pj( X 3L~ 



XjQ), 



3 



XjR-XjQ), 
3 



and 



(10) 



a q q 

the environment Hamiltonian can be rewritten as 

whereas the Hamiltonians H e and if p acquire the stochastic phase factors: 

H e = -A L a\ a Q - A R e* fl a ] R a Q + H.c., (12) 



H p = -A A e* A b\ b Q - A B e* B b ] B b Q + H.c. (13) 



3^1)1 
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£b = 'Y^Pj{X jB -X jQ ). 

j 

For simplicity, we assume that there are no phase shifts for the electron transitions between 
the electron source S and the site L, and the electron drain D and the site R, so that x jS = 
XjL, and XjD = XjR, with the same assumption for the N-A and P-B proton transitions, 
XjN = XjA and Xjp = Xj B . 

B. Rate equations 

The time evolution of the electron operators n a is determined by the Heisenberg equations: 

h L = iA L e <L a[a Q -i^2 c\s a L + H.c, 

k 

n R = iA R e <R a R a Q -i^2 f kD c{ D a R + H.c, (14) 

k 

and 

h Q = -iA L e* L a[ a Q - iA R e^ R a R a Q + H.c. (15) 
For the proton populations N a , we derive the similar set of Heisenberg equations, 

N A = %A A e* A b A b Q -tJ2 T qN d\ N b A + H.c, 

q 

N B = iA B b B b Q -iJ2 T qP d\ P b B + H.c. (16) 

i 

This set should be complemented by the equation for the proton population of the interaction 
site, 

N Q = -iA A e* A b\ b Q - iA B b B b Q + H.c, (17) 
as well as by the equations for the operators of electron and proton reservoirs, 

icks = £ks Cks — tks o<l, 

iCkD = £kD CkD — tkD dR, (18) 
% d q N = E q N d q N ~ T q N b A , 

id q p = Egpdgp -T qP b B . (19) 
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1. Contribution of reservoirs to the rate equations 

It follows from Eq. ( fl8|) that the electron operator Cks can be represented as 

c kS = cf s -t kS [ dtiH[cS(*).«S t (*i)]+>^(*i)«(*-*i). (2°) 



where c k °g(t) is the free variable of the S'-lead, and 9(t — t{) is the Heaviside step function. 
Similar expressions take place for the electron operator CkD{t) and for operators d q N, d q p of 
the proton reservoirs. For the weak coupling between the reservoir S and the electron site 
L we obtain 



(4(*)4°i(*)) = -ihs J dhicfP^cf^t)) ([a^a^t^eit-h). (21) 

Thus, contribution of the S'-lead to the evolution of the average electron population (n^) 
(see Eq. (TH|) ) is determined by the expression 

i £fo<«l(*)Q*s(f)) = "El^l 2 / Mi^c^tt)) (a{{t)a L { tl )) - 
k k J 

(^(tOcSWXaLCtOaiCt))}. (22) 

The correlator (c k g (ti)c k °s(t)) is proportional to the Fermi distribution function, fs(£ks) of 
electrons in the reservoir S, 

(ct ) s\t l )ct ) s (t)) = f s (e kS )e~^ ™, (23) 

where the Fermi function, 



fs(e) 



exp ( — ] + I 



is characterized by the electrochemical potential us and temperature T. We assume that the 
site L is weakly-coupled to the reservoir S and to the site Q, thus, we can use free-evolving 
operators, 

a L (t) = e- i ^ t - t ^a L (t l ), 
to calculate the corresponding correlation functions in Eq. (1221) . e.g., 

(a[(t)a L (t l )) = (n L (t))e^-^. 

Introducing the energy-independent rate constant, 

7s = 27iJ2\tks\ 2 S(e L -e kS ), (24) 



we calculate the contribution of the 5-lead to the time evolution of the population (til), 
i E*m (4(*)Pfcff(*)> + H - c - = 75 [/s(^) - (n L >]. (25) 

k 

The same analysis can be applied for a calculation of contributions of the electron lead D 
and the proton leads N and P to the corresponding populations (n R ) and (Na), (Np). The 
proton transfer rates between the sites A and C and the negative and positive sides of the 
membrane, respectively, are determined by the coefficients Yn and Yp where, e.g., 

Y N = 2tt J2 \ T in\ 2 S(E a - E qN ). (26) 
i 

2. Contribution of site-to-site tunneling to the rate equations 

To calculate a contribution of the L-Q tunneling to the evolution of the populations (rii) 
and (uq), we start with the amplitude ctQ, which obeys the equation 

ia Q = e Q a Q - u N Q a Q - A* L e~^ L a L - A* R e~^ R a R . (27) 

In the case of weak L-Q and R-Q tunnel couplings, the formal solution of Eq. (127]) can be 
written in the form 

I ^(-.^(^^^(tOl^lAle^^^a^tO + A^e^^^a^Xti)}, (28) 

where Oq (t) is the free operator of the site Q, obeying the equation (1271) with the tunneling 
terms neglected (Aj, = 0, A R = 0). 
Taking into account the formula, 

iA L (e^a[a^) = 

\A L \ 2 J dt 1 (ag )t (t 1 )ag ) (t)) ([e^(t)4(t), a L (h))^ 6(t - h), (29) 

which is similar to Eq. (j210 . we obtain 

iA L (e^ a{a Q ) = \A L \ 2 j dt^e**™ e^W) (a^x) a Q (t)) (a^h) 4(f)) - 

(e <faW e -W) {aQ{t) a U tl)) (4( t) flL(tl)>} . (go) 
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Dropping the label we assume that the time evolution of the operators ctQ in Eq. f[3"Uj) is 
calculated with the free-evolution formula, 

a Q (t) = e-^(^) aqih) - e -*«(*-'0 [i _ e -o(*-*i)] Nq ^ aQ ^y (31) 

For the free-evolving proton operator of the interaction site we obtain a similar expression 
b Q (t) = e -«*Q(t-ti) 6g ( fl ) _ e -^Q(^i) [1 _ e «o(t-ti)] ng(tl) 6q (^). (32 ) 

The influence of the environment on the electron tunneling between the sites L and Q, 
and between the sites R and Q, is determined by the correlators (er 1 ^^ e 1 ^®} and 
( e itUt) e -i€t(*i)^ where 

(e ^W e -W) = exp{-zAi(t - tx)} exp{-A L T(t - tx) 2 }. (33) 



The reorganization energy, X L , is defined as 



15| 



-^(xjL-Xjq) 2 . (34) 

i 

The electron reorganization energy A#, and the proton reorganization energies and A B , 
are defined in a similar way. In particular, 

,2 

2' 



5. Equations for populations of electron and proton-binding sites 

Consequently, we derive the system of rate equations for the average populations of the 
electron sites, 

n L + lsn L = Isfsi^L) + ®l, 
n R + 7£>n R = 7d/d(£h) + ®r, 

n Q = -$ L - $r, (36) 

and for the average populations of the proton-binding sites, 

n a + t n n a = t n f n {e a ) + $,4, 

N B + T P N B = T P Fp{E B ) + $ B , 

Nq = -$a - $b- (37) 
11 



Here 3> CT (er = L, R) and $ a (a = A, B) are the functions of the average electron and proton 
populations, respectively. In addition, due to a strong electron-proton Coulomb interaction 
on the site Q, the kinetic terms $ CT and $ Q depend on the correlation function, 

(K) = (n Q (t)N Q (t)) : (38) 
of the electron and proton populations on the site Q, 

$ CT = K a (e a -e Q + X a )(n Q )(l - n a ) - K a (e a - e q - A CT )(1 - n Q )(n a ) + 
{K c (e a - 6q + u + X a ) - K a (e a - sq + A CT )}(1 - n a )(K) - 
{K a (e a - E Q + u - \ a ) - K a (e a - E Q - \ a )}(n a }(N Q - K), (39) 

where K a (e) is the Marcus rate for electron transfer between the site a and the interaction 
site Q, 

*M = l^l ! ^5«p (~) ■ (40) 

The proton term $ Q is determined by the expression, similar to Eq. (13"9~|) . as 

$ Q = K a (E a -E Q + K)(N Q )(1 - N a ) - n a (E a -E Q - A q )(1 - N Q )(N a ) + 
{K a (E a -E Q + u + Ac) - Ka(£« -E Q + A q )}(1 - N a )(K) - 
{K a (E a -E Q + u - A a ) - K a (E a -E Q - A a )}(N a )(n Q - K), (41) 

where K a (E) is the proton Marcus rate for the transitions between the site a and the proton- 
binding site Q, 

*.m = l^l^exp (-J^) . (42) 
4- Equation for the electron-proton correlation function 

For the correlator, (K), of the electron (uq) and proton (Nq) populations of the interac- 
tion site, we derive the following equation 

(K) =F L + F R + F A + JF B , (43) 

where 

JF a = K a (E a -e Q + u - X a ){n a ){N Q - K) - 
Kai^o ~£q + uo + A (J )(1 - n a ){K), 
JF a = n a (E a - Eq + u - A a )(N a )(n Q - K) - 

K a (E a -E Q + u + A a )(l-N a )(K). (44) 
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5. Electron and proton currents 



Electron currents Is, Id and proton currents In, Ip are determined by an increase of the 
number of particles, electrons or protons, in the corresponding reservoir. In particular, a 
variation of the electron number in the drain lead gives a current 

Id = J t J2( c kD°kD) = lD[{n R ) - f D (e R )], (45) 

k 

whereas the proton current Ip is given by 

I ? = J t E<4pV> = r P [(N B ) - F P {E B )\. (46) 
i 

Here, 

7D = 2vr ^ \t kD \ 2 5(e R - e kD ), 
k 

and 

r P = 2tt J2\T q p\ 2 S(E B - E qP ) 

are the electron (7^) and proton (Tp) transfer rates between the electron site R and the lead 
D, and between the proton-binding site B and the P-side of the membrane, respectively. The 
multiplications of the particle currents introduced above by the electron or proton charges 
produce the standard electric currents. 



It follows from Eqs. ( I36II371) that, in the steady-state, we have the relations: 

{h a ) = , (N a ) = , 

so that 

$L + ®R = , 
$4 + ®B = , 

and 

Is = (d/dt) J2(c{ s c kS ) = -I D , 



I N = (d/dt) ^{d\ N d qN ) = -I P . 
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6. Quantum yield of the electron- driven proton pump 



The productivity of the proton pump is determined by a quantum yield, 



■ID 



(47) 



and by the power-conversion efficiency 77, 



rj = QY x 



(48) 



At the standard conditions, we have 



Hp — = V p + 60 meV = 210 meV, 



and 



Us — f^D = V e = 600 meV, 



therefore 



7] ~ 0.35 x gy. 



If a quantum yield Qy is of order one (or 100%), the power-conversion efficiency 77 may be 
as much as 0.35 (or 35%). 

C. Langevin equation 

For the redox-loop mechanism of a proton translocation through the membrane, the 
electron and proton sites, labelled by the letter Q, are attached to the shuttle: a molecule 
diffusing between the N and P sides of the membrane (see Fig. 2). This Brownian motion 
can be described by the one-dimensional overdamped Langevin equation for the coordinate 
x of the shuttle, 



We assume that the shuttle molecule moves along a line connecting the sites L and A, 
located at x — —xq, and the sites R and B, both having the coordinate x = xq. The borders 
of the membrane, at x — ±xq, are schematically shown in Fig. 2. In Eq. (I49p . ( is the drag 
coefficient of the shuttle, and £ is the Gaussian fluctuation force, which is characterized by 
the zero-mean value, (£) = 0, and the correlation function, 




(49) 



mm) 



2(T5(t - , 
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proportional to the temperature T of the environment. The diffusion coefficient D of the 
shuttle is also proportional to the temperature: D = T/(. The motion of the shuttle is 
restricted by the membrane walls, which are simulated by the confinement potential U c (x), 



U c (x) = U d0 <l 



+ 



exp 



X + x c 



(50) 



having the barrier height U c q, the width 2x c (x c > x$) and the steepness l c . 
The potential barrier U s (x), 



U s (x) = U, 



sO 



rxp ( - - Xs ) + :i 



x + x s 
< l xp I — j ) + 1 



(51) 



does not allow the shuttle with a non-zero charge q = Nq — rig (in units of |e|) to cross the 
lipid interior of the membrane. This barrier is determined by the height U s q, the steepness 



l s , and the width 2x s . 



III. RESULTS 



We solve the rate equations (I36II37P for the electron (n CT ) and proton (N a ) populations 
jointly with the equation ( 1431) for the electron-proton correlation function on the site Q, 
K = (uqNq). Our approach can describe two mechanisms of the redox-linked proton 
translocation across the membrane: (i) the static interaction site Q and (ii) the situation 
when the site Q diffuses between the sides of the membrane. The mechanism (i) roughly 
corresponds to the proton pump operating in cytochrome c oxidase (CcO) 3|-|6|, whereas the 



design (ii) can be attributed to the redox loop mechanism, w 



and proton transfers in the inner membrane of bacteria [7H12 1 . 



1. 


wi 


7- 
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rich is responsible for electron 



A. Static proton pump 

Here, we consider the mechanism (i), where the interaction site Q does not change its 
position (see Fig. 1). We assume that protons are transferred across the membrane, from the 
negatively charged side N, with an electrochemical potential hn, to the positively charged 
side P, having an electrochemical potential fip. All potentials and energies are measured in 
meV. 
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1. Parameters 



The difference of electrochemical potentials, A/z# = Hp — Hn, is determined by the 
following expression 

A/in = V P - 2.3 (RT/F) x ApH, (52) 

where V^, is the transmembrane voltage, R and F are the gas and Faraday constants, re- 
spectively, T is the temperature (in Kelvins, ks = 1), and the concentration gradient ApH 
is about — 1. 

ii a. The coefficient 2.3 (RT/F) is about 60 meV at room temperature, 
T = To = 298 K. It follows from Eq. (!52l) that the potentials of the N and P sides of the 
membrane can be written as 

»n = -»m - AV 2 " 30 x (AT /T ), 

Hp = Hm + AV p /2 + 30 x (AT/T ), (53) 

where AV P — V p — Vq, AT = T — T . At the standard conditions, when T = T , V p = Vq = 
150 meV, for the electrochemical potential fino we have: fino — 105 meV. Thus, the total 
oroton gradient across the membrane, A/z#, is about 210 meV. As in the CcO proton pump 
3|, la], we assume that the proton-binding sites A, Q p , and B are located approximately on 
the line connecting the N and P sides of the membrane with the following coordinates: 
%a = 0.1, xq = 0.3, xb = 0.5. The coordinates of the sites are counted from the middle of 
the membrane in a direction towards the P-side and are measured in units of the membrane 
width W with W ~ 4 nm. Protons are delivered from the iV-side to the site A by the 
so-called D-pathway crossing about a half of the membrane. We also note that the B-site 
is located next to the P-side (see Fig. 1). An influence of the transmembrane voltage V p on 
the energy levels of the proton sites is described by the formulas 

E A = E AQ + x A x AV, 
E Q = E Q0 + x Q x AV, 

E B = E B0 + x b x AV. (54) 

For the proton energy levels, E Ao ,Eq , and E B o, at the voltage V p = V , we assume the 
following values (in meV): E A o = — 155, Pqo = 250, and Ebo = 185, unless otherwise 
specified. This means that at the standard conditions, the proton begins its journey at 
the iV-side with the potential /ijy = —105 meV and jumps to the A-site having a lower 
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energy (—155 meV). However, the next proton-binding site Q p has a much higher energy 
(~ 250 meV), so that the proton transfer cannot occur without a mediation of the electron 
component. The electron site Q e is electrostatically coupled to the proton-binding site Q p 
with the Coulomb energy uq. Thus, in the presence of an electron on the site Q e the energy 
of the Q-proton decreases to the level Eq — u ~ —220 meV, provided that u ~ 470 meV. 
Now the proton can move from site A to site Q, since Eao > Eqq — uq. Depopulation of the 
electron site Q returns the energy level of the Q-proton to its original value Eq$ = 250 meV, 
which is higher than the energy level of the next-in-line B site, Ebo = 185 meV, and is much 
higher than the energy level of the A-site. We assume that the backward proton transfer 
(from Q p to A site) is described by the inverted region of the Marcus formula, so that the 
probability of such transfer is low, compared to the probability of the proton transfer from 
the site Q p to the site B. No additional gate mechanism is necessary here. 

For the sake of simplicity, we assume that three electron-binding sites L, Q e , R as well 
as the source and drain leads are positioned on a line, which is parallel to the surface of 
the membrane (see Fig. 1). Thus, the transmembrane gradient V p has no effect on electron 
transport from the electron source S to the drain D. For the potentials of the electron 
reservoirs, we choose the following form 

/AS = /^eO + V e /2, 

= HeO ~ V e /2, (55) 

with yU e o = —500 meV and with the electron voltage gradient V e = 600 meV, unless otherwise 
indicated. The electron voltage gradient V e roughly corresponds to the drop of the redox 
potential along the electron transfer chain in the cytochrome c oxidase [i s)]. We assume 
that the electron pathway includes the source reservoir (/is = —200 meV), the site L (sl = 
—210 meV), the interaction site Q e (eq = —250 meV), the site R (er = —770 meV), and 
the electron drain reservoir having the potential /ip = —800 meV. 

We assume that the electron and proton transfer between the active sites, L-Q, R-Q and 
A-Q, B-Q, are quite fast, with amplitudes A L ~ A# ~ 0.3/ps and ~ A B ~ 0.3/ps, 
whereas the transitions to and out the electron and proton reservoirs are characterized 
by much slower rates: 75 ~ 7^ ~ 1.5/ns, and — — 0.75/ns. The responses of 
the environment to the electron and proton transitions are described by the corresponding 
reorganization energies: = \r = X e and = A B = A p , respectively. Here, for the 
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standard case, we assume that A e ~ 100 meV and A p ~ 100 meV. This set of parameters 
provides an efficient operation of the redox-linked proton pump. 

2. Dependence of the proton current on the transmembrane voltage 

In Fig. 3, we show the steady-state proton current Ip as a function of the transmembrane 
voltage gradient V p , at three different values of the electron voltage: V e = 500, 600, 700 meV. 
We use here the standard set of other parameters (see the previous subsection), where 
T = 298 K and A e = A p = 100 meV. 

The proton current Ip is equal to the number of protons pumped energetically uphill 
(at V p > 0), from the negative side N to the positive side P of the membrane, per one 
microsecond. At the difference V e = 600 meV of source and drain redox potentials, the 
system pumps more than 200 protons per one microsecond against the transmembrane 
voltage gradient V p = 150 meV. According to Eq. ( 1531) . this voltage corresponds to the proton 
electrochemical gradient A/i# = 210 meV, which is usually applied to the internal membrane 
of mitochondria and the plasma membranes of bacteria. The number of pumped protons goes 
down as the proton voltage V p increases, and goes up with increasing the electron voltage 
difference V e . The proton current saturates at V e > 750 meV. It is evident from Fig. 3 
that at high enough electron voltages (V e > 600 meV), the pump is able to translocate 
more than 100 protons per microsecond against the proton gradient V p , exceeding 250 meV 
(A/i > 310 meV). The quantum yield QY is about one (with a power-conversion efficiency 
r\ ~ 35%) in the whole region of electron and proton voltages: 500 meV < V e < 800 meV, 
< V p < 300 meV. 

3. Proton current and the quantum yield as functions of temperature 

Figure 4 shows the pumping proton current, Ip (i.e., the number of protons translocated 
from the negative to the positive side of the membrane per one microsecond) versus the 
temperature T measured in Kelvins. The graphs are presented at three values of the electron 
and proton reorganization energy: A = 100, 150, 200 meV. We assume here that A e = A p = 
A, with the electron voltage V e = 600 meV and the proton gradient V p = 150 meV. It is of 
interest that at A > 150 meV the pumping current has a pronounced maximum near the room 
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temperature, 200 K < T < 300 K, although the quantum yield is higher, QY ~ 1, at lower 
temperatures. The performance of the pump deteriorates at higher reorganization energies 
when the coupling to the environment increases. Increasing the reorganization energy leads 
to increasing the probability for an electron to be transferred through the system, losing 
all its excess energy to the environment without transferring this energy to protons. Such 
a probability is further increased at large temperatures leading to the observed decrease of 
the quantum yield. 

4- Dependence of the proton current on the parameters of the interaction site 

The energy transfer from the electron to the proton component occurs on the interaction 
site Q = {Q e ,Qp}, which has one electron (eq) and one proton (Eq) energy levels (see 
Eq. f)54p ). The electron on the site Q e is electrostatically coupled to the proton, which 
populates the site Q p , with the Coulomb energy w - It follows from Fig. 5 that the proton 
pumping current Ip exhibits a resonant behavior as a function of the charging energy uq and 
the position of the proton energy level Eqq. The dependence of the pumping current on the 
electron energy eq has a resonant character as well. Here we assume that V e = 600 meV, 
V p = 150 meV, A e = A p = 100 meV, and T = 298 K. The energetically-uphill proton current 
has a pronounced maximum (Ip ~ 220//is) at the Coulomb energy w = 470 meV and the 
proton energy Eq = 250 meV, provided that the electron energy Eq = —250 meV. It is 
important that the proton pump is robust to the variations of the Coulomb energy uq and 
the proton energy Eqq in the range ±50 meV from the resonant values. The quantum yield 
QY is very close to one in the central region of Fig. 5, so that the power-conversion efficiency 
i] is about of 35%. 

Figures 3, 4 and 5 clearly demonstrate that, at standard physiological conditions, the 
static redox-linked proton pump ( "CcO-pump" ) efficiently converts the energy of electrons to 
the more stable energetic form of the proton electrochemical gradient across the membrane. 

B. Redox loop mechanism of electron and proton translocation 

In many biological systems, electrons and protons can be transferred across a membrane 
by means of a molecular shuttle diffusing inside of the membrane, from one side to another. 
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Here we show that the mathematical model described in Section II can be successfully 
applied for a description of the redox loop mechanism, which utilizes the Brownian motion 
of the shuttle Q carrying both electron, Q e , and proton, Q p , sites (see Fig. 2). As in the 
previous case, we have to solve here a system of master equations for the electron (til, tiq, tir) 
and proton (Na, Nq, N b ) populations, Eqs. (I36|37p . and for the correlation function K of 
electron and proton populations on the site Q, Eq. ( H3|) . However, these master equations 
should be complemented by the Langevin equation, Eq. fj49|) . for the time-dependent shuttle 
position x. We note that the electron tunneling between the sites L-Q, Q-R, as well as 
the proton transfer rates between the sites A-Q and Q-B, depend on the position x of the 
shuttle. 

1 . Parameters 

We assume that the electron site L is located near the negative (N) side of the membrane, 
at where Xq = 2 nm. The other electron site R is near the P-side of the membrane, 

at x = +Xq. The reservoir S, connected to the site L, serves as a source of electrons, 
and the reservoir D, coupled to the site R, serves as an electron drain (see Fig. 2). The 
tunneling amplitudes A L ,A R are determined by the amplitudes A L0 , A R0 , and by the 
electron tunneling length l e : 



The proton-binding site A is located at the end of the iV-side proton pathway, whereas 
the site B terminates a pathway, which goes into the P-side of the membrane. For the 
x- dependencies of the proton transfer amplitudes and A^, we choose the following 
relations: 



where l p is the proton transfer length. It should be noted that our model produces the 
same results when the proton amplitudes are given by the expressions similar to Eqs. fl56|) . 




(56) 




(57) 
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For the transfer parameters, we choose the following values: A L0 ~ A R0 = 0.04 meV, 
A,4o ~ A so = 0.04 meV, and l e = 0.25 nm, l p = 0.25 nm. Couplings to the electron and 
proton reservoirs are described by the rates 75 ~ 7d = 0.5/ns and T N ~ T P = 0.1/ns. The 
system is robust to significant variations of the transfer parameters. 

The confinement potential U c (x) is determined by the height U c = 500 meV, the steepness 
l c = 0.1 nm, and the half-width x c = 2.7 nm. The potential barrier U s (x), preventing 
the charged shuttle from entering into the membrane, is characterized by the height U s = 
770 meV, the width x s = 1.7 nm, and the steepness l s = 0.05 nm. 

Accordingly, the electron and proton populations of the shuttle are almost completely 
compensated, hq ~ Nq, so that the potential U s (x) gives a negligible contribution to the 
energies of electrons and protons. However, we have to take into account the fact that in the 
presence of the voltage gradient, V p ~ 150 meV, the electron (eq) and proton (Eq) energies 
on the moving shuttle depend on the shuttle position x: 

X M 

E Q = E Q0 + ^-V p , (58) 

ZXq 

with Eq = 280 meV, and Eq = u /2 = 200 meV, where for the charging energy u of the 
shuttle we have: u = 400 meV. 

Thus, electrons move from the source reservoir, having the electrochemical potential 
fis = 420 meV, to the L— site (with the energy e L = 380 meV), and, thereafter, to the 
shuttle. On the opposite side of the membrane, the electron, populating the shuttle, jumps 
to the site R (er = —170 meV) and, finally, to the drain reservoir ([ip = —230 eV). 
The total drop of the redox potential in this electron-transport chain can be estimated as 
fJ'S — = 650 meV. 

Protons move from the iV-side of the membrane (fx N = —105 meV) to the site A, having 
a lower energy Ea = —150 meV. The energy level Eq = 125 meV of the proton on the 
shuttle, located near the iV-side of the membrane (at x = —xq), is much higher than Ea, if 
the shuttle contains no electrons. However, the shuttle populated with a single electron is 
more attractive for protons, since in this case the effective energy of the proton, Eq — u = 
—275 meV, is less than the energy of the proton-binding site A. The shuttle, carrying one 
electron and one proton, diffuses to the opposite side of the membrane (x = +x ), where the 
electron, with energy eq — u = —195 meV, is able to tunnel to the site R, having a slightly 
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higher energy er = —170 meV. In the absence of an electron, the energy of the proton on the 
shuttle (at x = +x ) increases to the level Eq = E Q0 + V p /2 = 275 meV, which exceeds the 
energy of the proton on the site B: E B = 150 meV. Consequently, the proton moves from 
the shuttle to the site B and, thereafter, to the P-side of the membrane characterized by the 
electrochemical potential fip = +105 meV. Thus, this redox loop mechanism translocates 
protons across the membrane against the proton electrochemical gradient A/ih = Hp — ^n = 
210 meV, and against the transmembrane potential V p ~ 150 meV. 

2. Proton translocation process 

Figure 6 exhibits the electron and proton populations of the shuttle, riQ(t) and Nq(t), 
correlated with the shuttle's position x(t) at T = 298 K, V p — 150 meV, and at A/x = 
210 meV. In this figure, we also show the time dependencies of the number of electrons, 
no(t), transferred to the drain reservoir, and the number of protons, Np(t), translocated to 
the positive side of the membrane. The shuttle diffuses between the membrane walls located 
at x — ±x (xo = 2 nm) with an average crossing time At ~ 2.5 /is. This time-scale is 
closely related to the diffusion time, 

t D ~ (Ax 2 )/2D ~ 2.66 fjs, 

obtained at a/ (Ax 2 ) ~ 2x = 4 nm, for the diffusion coefficient of the quinone molecule 
D ~ 3- 10~ 12 m 2 /sec. 

At t ~ 0, the shuttle, located at x ~ — xo, is loaded with one electron and one proton 
taken from the negative side of the membrane (see Fig. 2). When t ~ 2.5 /is, the shuttle 
reaches the positive side (x = +xq = 2 nm) and unloads the electron to the the site R 
(and later to the drain lead D) and the proton to the site B, coupled to the P-side of 
the membrane. Consequently, the population N P of the P-side grows. The empty shuttle 
diffuses back, to the iV-side, completing the cycle, and the process starts again. In twenty 
microseconds, the shuttle performs four complete trips and translocates about four electrons 
and four protons across the membrane. 
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3. Voltage and temperature dependencies 

The numbers of electrons and protons, rir> and Np, respectively, transferred across the 
membrane in one millisecond, are shown in Fig. 7 as functions of the transmembrane proton 
voltage V p . The electrochemical gradient of protons, A/i = /j, P — fi N , is proportional to 
V p : A/i ~ V p + 60 meV (at T = 298 K). The results in Fig. 7 are averaged over ten 
realizations. The system is able to translocate more than 120 protons per ms against the high 
transmembrane voltage, V p < 250 meV, that corresponds to the electrochemical gradient 
Afj, < 310 meV. 

It follows from Fig. 8 that the translocation mechanism works efficiently in a wide range of 
temperatures, 250 K < T < 500 K. In this range, the system pumps more than 120 protons 
per millisecond with a quantum yield exceeding 90% and with a power-conversion efficiency r\ 
higher than 40%. With increasing temperature, the shuttle performs more trips between the 
sides of the membrane, thus, carrying more electrons and protons. This increases the proton 
current (i.e., the number of protons translocated per unit time). We note that the proton 
population of the shuttle occurs only after loading the shuttle with an electron. At very 
high temperatures, T > 500 K, the shuttle moves quite fast, and protons have less chances 
to jump on the shuttle. Consequently, the gap between electron and proton currents grows 
with the temperature, thus deteriorating the performance of the pump. 

IV. CONCLUSION 

Two different mechanisms of energetically-uphill proton translocation across a biomem- 
brane are described by the same physical model. This model includes three redox sites 
(L, Q e , R) and three proton binding sites (A, Q p , B) attached to the source (S) and drain 
(D) electron reservoirs, as well as to the proton reservoirs on the positive and negative sides 
of the membrane. We have shown that it is the strong Coulomb interaction between the 
electron site Q e and the proton site Q p , which plays the most prominent role in the process 
of energy transformation from electrons to protons. In this case, the whole electron-proton 
transport chain can be divided into weakly coupled clusters of sites, so that the total num- 
ber of occupation states is equal to the sum (not to the product) of occupation states in 
each cluster. At physiological conditions, our model demonstrates a proton pumping effect 
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with a quantum yield near 100% and a power-conversion efficiency of order of 35%, for both 
the static proton pump, related to the cytochrome c oxidase, as well as for the redox-loop 
mechanism, where electrons and protons are translocated by the diffusing molecular shuttle. 
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FIG. 1: (Color online) Schematic diagram of the static proton pump. The electron transport chain 
starts at the source (S) lead. Thereafter, high-energy electrons, e~, tunnel energetically- downhill 
(through the yellow path) to the sites L, Q e , R and, finally, to the drain D. Low-energy protons, 
H + , move energetically-uphill (in blue) from the negative (N) side of the membrane to the sites 
A, Q p , B and, eventually, reach the positive (P) side of the membrane. 
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FIG. 2: (Color online) Schematic diagram of the redox loop mechanism. Here, the electron- 
proton interaction site, Q = {Q e ,Qp}, is placed on the molecular shuttle (shown in green), which 
diffuses along the line connecting the negative and positive sides of the membrane. From the source 
reservoir S, an electron e~ jumps to the site L and, thereafter, to the shuttle, located at x = —xq. 
The shuttle also accepts a proton H + transferred from the iV-side of the membrane via the site 
A. The loaded shuttle moves randomly toward the positive side (P) of the membrane, where (at 
x = xq) the electron is subsequently transferred from the site Q e to the site R and to the drain 
reservoir D, and the proton jumps from the site Q p to the site B and, finally, to the positive (P) 
side of the membrane. We note that, in this design, the electron site L and the proton site A 
are located near the iV-side of the membrane (shown by the horizontal blue dashed line) , and the 
electron site R and the proton site B are placed near the P-side. 
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FIG. 3: (Color online) Proton current versus transmembrane voltage V p at room temperature, 
T = 298 K, and three different electron potentials: V e = 500, 600, and 700 meV. The proton 
current is almost constant for low values of V p , and decreases for increasing V p . 
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FIG. 4: (Color online) Proton current (the number of protons translocated across the membrane 
per one microsecond) and quantum yield versus temperature for the electron voltage V e = 600 meV, 
transmembrane proton voltage V p = 150 meV, and three different reorganization energies: A = 
100, 150, and 200 meV. The proton current and quantum yield both decrease, for increasing A. 
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FIG. 5: (Color online) Dependence of the proton current (the number of protons pumped across 
the membrane per one /xs, see the color bar on the right side) on the charging energy uq, and on 
the energy Eqq of the central proton site for V e = 600 meV, V p = 150 meV, and T = 298 K. 
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FIG. 6: (Color online) Time evolution of the electron-proton translocation process. Here x is 
the location of the shuttle, tiq and Nq are the electron and proton populations of the shuttle, 
respectively, rip is the number of electrons transferred from the electron source S to the electron 
drain D, and Np is the number of protons translocated from the negative (iV) to the positive (P) 
side of the membrane. It can be seen from this figure that the loading/unloading of the shuttle 
with electrons and protons, as well as the electron and proton transfer across the membrane, are 
clearly correlated with the spatial motion of the shuttle. 
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FIG. 7: (Color online) Numbers of electrons, np, and protons, Np, translocated across the 
membrane in one millisecond, versus the transmembrane proton voltage V p at room temperature, 
T = 298 K, and at (yu<j — fio) = 650 meV. Clearly, it is much harder to transfer protons against 
the higher transmembrane voltages. 
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FIG. 8: (Color online) Temperature dependence of the numbers of electrons, rir>, and protons, 
Np, transferred across the membrane by the diffusing shuttle, at V p = 150 meV and (fj,g — fio) = 
650 meV. We also present here the quantum yield, QY, and the power-conversion efficiency, rj, 
of the process as functions of the temperature. At higher temperatures, the shuttle moves faster 
and carries more electrons and more protons. However, if the temperature is too high, the shuttle 
has not enough time to be loaded with electrons and protons, and sometimes travels empty. As a 
result of this, the electron and proton currents decrease at high temperatures, thus decreasing the 
efficiency of the pump. 
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